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SUMMARY 


Based on photogeologic mapping of Viking Orbiter images of Mars, four 
distinct informal stratigraphic units can be defined for the region north 
of 70°N latitude. They are: (a) bulbous plains, (b) mantled plains, (c) dune 
deposits, and (d) layered deposits/perennial ice. 

The bulbous plains unit underlies all other north polar units and 
represents a sub-unit of the mottled cratered plains. Based on crater size- 
frequency data, bulbous plains is equivalent in age to the relatively 
old cratered plains unit at Tempe Plateau (90°W, 30°N). The low albedo of 
bulbous plains and the appearance of what appears to be a dike suggest that 
bulbous plains has a volcanic origin. Cumulative crater size- frequency 
distribution functions for bulbous and mantled plains display two-segment 
curves with a crater-production slope of -2.0 and a cratering-cbl iteration 
equilibrium slope of -0.7. This supports the interpretation that mantled 
plains were formed by dust blanketing of bulbous plains. As calculated from 
crater size-frequency distribution functions, the relative surface-obliteration 
rates are 1:9:88 for bulbous, moderately mantled, and heavily mantled plains, 
respectively. 

The layered deposits probably represent a facies of mantled plains 
in which dust has been deposited onto the perennial ice cap rather than directly 
onto the ground. Thus, the areal extent of the layered deposits at the north 
and south poles may indicate the maximum extent of the perennial ice caps as 
controlled by changes in Martian orbital parameters. 

Dune deposits occur in the form of longitudinal, transverse, and barchan 
dunes, and possibly as sheet sand deposits. Actual dunes are generally confined 
to mantled plains, suggesting that mantling provides a proper substratum for 
dune accumulation. 

Surface winds for the north polar region were determined from dune 
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orientations to flow predominantly counterclockwise around the polar cap, 
although there Is evidence for dune modification by secondary winds spiralling 
clockwise off of the cap. 

A gradation from densely spaced transverse dunes to more dispersed 
barchan and transverse dunes implies a general thinning of circumpolar dune 
deposits downwind of extensive areas mapped as bulbous plains. The author 
suggests that dune material is being stripped from bulbous plains and Is 
accumulating as dunes downwind on mantled plains. 
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Telescopic Observations 

Man first became interested in the polar regions of Mars when white 
polar spots were reported by Cassini in 1666 and again by Huygens in 1672. 

In 1784, Sir William Herschel studied the changing size of the polar caps 
with seasons and interpreted the polar spots as seasonal accumulations of 
ice and snow. The melting of the polar caps was found to be intimately 
related to the "wave of darkening" of the Martian surface during summer (1). 

Although Wilhelm Beer and Johann Madler reported the presence of a dark 
band around the receding north pole of Mars in 1830, little further attention 
was given to the observation until 1894 when Percival Lowell reported a 
similar band around the south polar cap. This sparked controversy as to 
whether the dark bands were real features of the Martian surface or just 
an optical illusion (2). Since the circumpolar bands were interpreted as 
moistened ground, their presence would imply actual melting of the retreating 
polar water ice rather than sublimation (1). Lowell envisioned his Martian 
"canals" as an irrigation system designed to bring water from the melting 
caps to the arid lower latitudes (2). 

Two distinct types of Martian polar "caps" have been recognized since 
1907. The spring and summer caps are surface deposits whereas the "caps" 
observed in autumn and winter are cloud mantles (i.e. the polar hood). The 
polar hood begins forming as the summer ends and disappears rather abruptly 
at vernal equinox (1). 

Two major differences between the north and south polar ice caps have 
been recognized from telescopic observations. Whereas the north polar cap 
is centered almost precisely on the rotation axis of Mars, the center of the 
south polar cap deviates from the rotation axis by 5°. Also, the north polar 
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cap is larger than the southern cap, probably due to the fact that the subsolar 
latitude during perihelion is presently in the southern hemisphere. It is 
probable that during this period of excess insolation in the southern hemi- 
sphere, the north polar cap grows at the expense of the southern cap (3). 

Mariner Mission Results 

Neither polar cap was imaged by the 20 pictures returned by the Mariner 4 
Mars flyby mission in 1965. However, spectrometry data taken during Mars 
occultatiori determined that the Martian atmosphere was essentially carbon 
dioxide at <a few millibars. This led many to argue that the polar caps may 
also be primarily carbon dioxide (2). 

The first "close-up" images of the south polar cap were taken by the 
Mariner 7 spacecraft as it flew by Mars in 1969. No dark circumpolar band 
was otserved for the south polar region, placing serious doubt on the possibility 
that a dark "melt band" might occur around the north pole (2). 

The Mariner 9 orbiting spacecraft photographed Mars from November 1971 
to March 1972, and provided excellent south polar coverage. This was, in 
part, due to the fact that the south pole was the only feature clearly seen 
through the global dust storm which enshrouded Mars during the early mission, 
and because the orbit of the spacecraft traversed the south pole every rev- 
olution during the orbit's periapsis. The entire south pole was imaged at 
3-4 km resolution with 20% covered by 300-400 m resolution (3). 

One of the most striking features of the south pole is the "quasi-! inear, 
swirl-like" pattern seen in the ice cap. The north polar cap was later 
found to display an even more pronounced swirl pattern (see figure 1). Based 
on Mariner 9 images, several mechanisms have been proposed to account for these 
features. Murray and Malin (4) proposed that the curvilinear features repre- 
sented the edges of stacked plates of ice and suggested that an offsetting 
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of these plates provided evidence for polar wandering. Others believed the 
features represented troughs rather than the edges of plates, and suggested 
spiralling polar winds (5), or mass-wasting (6) as the dominant erosion 
mechanisms. 

The excellent Mariner 9 coverage of the south polar region provided a 
clear understanding of its stratigraphic sequence (7). The ancient cratered 
terrain, which is the oldest unit on Mars and covers most of the southern 
iirnvi cohere, "^orl ies the south polar region as well. Massive (i.e. un- 

blankets, probably composed of fine, wind-transported dust, overlie 
tee -■'"•ciont cratered terrain in the south polar region. These blankets have 

Pensively eroded to form pitted terrain and etched terrain (3). The 
massive blankets are overlain, sometimes unconformahl y, by layered deposits 
which Ivr'o been subsequently eroded to expose the laye“i^g in troughs and 
scarps. Tn addition to the thin mantles of C0 ? ice which advance o-vj ^etroa*. 
each year (7), permanent ice caps cover both poles though Ma’-Wnr.- o inves- 
tigations could not determine whether the permanent caps consist r>r iyi <v* 

CO? ice. 

The layered deposits probably form by cyclic deposition of wind-blown 
dust carried in from the mid-latitudes (9). Deposition models using 50,000 
and 95,000 year climatic cycles have been proposed to account for the layers 
(4), although Cutts concludes that such a period is too short to allow transport 
of sufficient material to the poles. However, a 2 million-year cycle of 
orbital ellipticity of Mars may modify the intensity of global dust storms 
enough to cause the layering, and allow time for a sufficient amount of 
material to be transported to the poles (9). 

Compared to the south polar coverage. Mariner 9 coverage of the north 
pole of Mars was relatively poor. Not only was the spacecraft's orbit much 
higher over the north pole compared to the orbital range during south polar 
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traverses, but the polar hood obscured the north polar region during most of 
the standard mission. Although visibility improved greatly during the Mariner 9 
extended mission, north polar coverage was severely limited because of a shortage 
in the spacecraft's attitude-control gas. 

Mariner 9 imagery products showed the north polar cap and its associated 
layered deposits to be very similar to deposits of the south pole. Minor 
differences include the more pronounced "swirl-like" pattern of the north 
polar cap and the greater extent of the northern permanent ice cap. Layered 
deposits in the north are more extensively covered by permanent ice than at 
the south pole. 

Although the lack of a dark circumpolar band in the south had placed 
serious doubt on the presence of such a band in the north, a very low albedo 
band around the northern pole was confirmed by Mariner 9. However, the feature 
was found to be fixed in position, rather than a halo surrounding a retreating 
ice cap. Sagan et. al. (10) attributed the presence of the dark band to the 

t 

scouring of a thin layer of bright dust by strong polar winds. 

I 

Mariner 9, for the first time, allowed determination of the stratigraphic 
sequence of the north polar region (11). Like the south polar region, the 
sequence in the north was found to include an underlying cratered unit, a i 

massive deposit interpreted as wind-blown dust, and layered deposits capped 
by permanent ice. 

However, the underlying moderately cratered unit is not the ancient 
cratered terrain, but rather a younger unit designated as mottled cratered 
plains. These plains are characterized by Soderblom et. al. (11) as having 
an abundance of small craters (10-20 km diameter), lower average albedo than 
most of Mars, and a highly mottled appearance which results from numerous 
bright and dark streaks and bright material highlighting crater rims and 
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interiors. As mapped by Soderblom et. al. the mottled cratered plains are 
distributed rather extensively between 55° to 80° north latitude*. 

Overlying the mottled cratered plains is a bright debris mantle which 
was subdivided into smooth, etch-pitted, and rippled plains (11). Rippled 
and etched-pitted plains were interpreted as depositional end erosional 
modifications of smooth plains. Because of the severely limited extent and 
poor quality of the Mariner 9 north polar coverage, the geologic mapping of 
Soderblom et. al . (11) is labelled by the authors as "preliminary" and contains 
numerous implied boundaries between units. 

V i In ng Mi ssio n Results ( Augu st. 1976-November 1978 ) 

Compared to Mariner 9, the Viking Orbiter spacecrafts, particularly 
V02, provided much improved coverage of the north polar region of Mars. 
Following the successful landings of Viking l.anders 1 and 2, the inclination 
of the orbit o' Viking Orbiter 2 was changed from 55° to 71° in nrdo~ to 
permit extensive north polar coverage before polar hood formation (33). 

Ground resolution of the north polar images taken with the VO? in this o**hit 
wa.> 100-180 meters as compared to 300-400 motors for Marine" 0 high-resolution 
coverage of the south pole (3). During the extend Viking mission, the orbital 
inclination of V02 was further increased to 00" and the periansis lowered 
from 1500 km to 778 km (33). This improved the ground resolution of the north 
polar imagery to about 60 meters. 

In addition to the Viking Orbiter imaging capabilities, the Orhiters 
included an infrared thermal mapper (IRTM) which provided conclusive 
evidence that the residual ice cap at the north pole was composed of a "dirty 
water ice" (12). Summer temperatures at. the south pole have not dropped low 
enough to allow a similar determination of the composition of the south polar 
cap (34). 
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The spiral patterns within the polar caps have continued to intrigue 
Viking investigators. Although it. is now generally accepted that these features 
are troughs rather than plate edges, there is still much controversy as to the 
origin of the troughs. As an alternate hypothesis to the rotating wind 
erosion and mass wasting processes mentioned earlier, Howard (13) has suggested 
that isolation-controlled ice ablation and accumulation could account for the 
spiral pattern of the troughs, with wind action only serving to redistribute 
dust across layered deposits. 

The circumpolar dark band, which had been confirmed by Mariner 9 
investigations, was determined from Viking images not to be a deflation 
feature as had been predicted earlier, but corresponded instead to large dune 
fields partially surrounding the north polar ice cap (14). It will be shown 
in this paper, however, that not all of the dark areas can be explained by 
the presence of visible dunes. 

Determining the source of such a large supply of dune material has 
become a very intriguing problem. Any model attempting to define the source 
of the north polar dune deposits must consider: (a) the unique geographic 
association of large fields with the north polar region, (b) the very low 
albedo of the dune deposits (A^ <_ 0.20) (12), and (c) the particle size 
required for saltation in the Martian environment (diameter = 200-250 um; 15). 

An added complication to the grain size requirement is that sands appear to 
be absent at the Viking Landing sites, possibly due to self-destruction of 
saltating grains during the impacts with rocks (i.e. the "kamikaze" effect; 

16). Therefore, the model mjst account for sand-sized particles being 
present at high latitudes while apparently absent at lower latitudes. Source 
models proposed fall into three general categories: (a) dune sands have been 
transported from lower latitudes predominantly by saltation (14), (b) dune 
material consists of aggregates of atmospherically-derived dust, which have 



9 

been deposited either directly onto the ground as "dirty iceballs" (17), or onto 
the polar cap where it is subsequently eroded from the layered deposits (14), 
and (c) dune sands have a regional source implying processes unique to the 
north polar region (15, this paper). Each model proposed has certain problems 
which must be accounted for. 

Blasius et. al. (35) have suggested that aeolian carving of many of the 
channels leading into Chryse Planitia may have supplied both suspended and 
saltating material to the northern plains. Although conditions may have 
been different in the past, neither wind streak data (J, 10) nor atmospheric 
models based on present corditions (18, 19) support such an atmospheric 
flow direction. 

A major complication with the proposed atmospheric-dust models is that 
they require aggregation of fines to form sand-sized particles, using 
aggregate honding mechanisms such as electrostatics (20), ice (17), or 
evanorit.es (14). A more serious problem is that the atmospheric-dust models 
fail to account for the low albedo of the dune deposits, since atmospheric 
dust on Mars is commonly known to have a high alhodo. 

A regional -source model must explain the uniqueness of the large supply { 

of dune material in the north polar rpgion, by assiming either a unique J 

source unit, a unique polar process, or both. Sguvros (15) has suggested that * 

a regional source of "coarse-grained" basalt may produce low-albedo dune 
material which retains its sand-size and low albedo because of low tJV weathering 
rates. Retention of an adequate sand supply may be aided by high mechanical 
weathering rates which result from freeze-thaw processes and glacial grinding. 

Based on the improved coverage by the Viking Orbiter spacecrafts, Squyres 
(15) revised the stratigraphic sequence in the north polar region to include 
dune material. Following much the same sequence as set. forth by Soderblom et. al.. 
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Squyres divided the stratigraphy into moderately-cratered plains material, 
debris mantle (massive deposits), dune material, laminated (layered) 
deposits, and perennial ice. Using results derived through photogeologic 
mapping and from crater size-frequency distribution data, it will be shown 
in this paper that the stratigraphic mapping of Squyres is inconsistent 
with his description of the units. 

The results of the present paper were primarily obtained from photo- 
geologic mapping of Viking Orbiter imagery in the form of mosaics and individual 
enlarged frames. The goals of the research were (a) to redefine and map 
distinct informal stratigraphic units in the north polar region of Mars 
(i.e. > 70°N latitude), (b) to define and interpret crater size-frequency 
distribution functions for particular units, (c) to map the directions of 
dune-driving winds based on the orientation and form of circumpolar dunes, 
and (d) to infer a source for the north polar dune deposits based on the 
relationships between dune-driving winds and the distribution of stratigraphic 
units. 
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Figure 1: Mariner 9 mosaic of the north polar region of Mars (from Mars 
as Viewed from Mariner 9 , NASA SP-329, 1976). Note the well developed 
spiral pattern within" the ice cap. The low-albedo band surrounding the 
cap corresponds in part to the presence of dune fields. 
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Flaure 2: Stratigraphic map of the north polar region of Mars: (1) dunes 

(2) layered deposits/perennial ice, (3) mantled plains, and (4) bulbous 
plains. Mantled plains include moderately- and heavily-mantled plains; 
bulbous plains Include bulbous and lightly-mantled bulbous plains. 


— „ 






DESCRIPTION OF STRATIGRAPHIC UNITS 


Bulbous Plains 

The most heavily-cratered terrain In the north polar region is informally 
referred to as bulbous plains by this author. It is characterized by an 
abundance of, small craters (<2 km diameter) relative to other north polar 
units (figure 9)» by numerous debris-flow craters, and by a mottled, bulbous 
appearance (figure 3). The mottled, bulbous appearance is predominantly a 
result of bright material surrounding dark, highly irregular polygons, which 
are 3-8 km wide. Mottling is also due to concentrations of bright material 
inside and on the exterior rims of many craters. Within the same Viking 
Orbiter frame, the albedo of bulbous plains is generally intermediate 
between very dark dune deposits and bright mantled plains. 

Mantled Plains 

Mantled plains display a higher albedo than either bulbous plains or 
dune deposits, and has a crater population characterized b.v a deficiency 
of small craters relative to bulbous plains (figure 9) and by shallow craters 
cornnonly lacking visible ejecta (i.e. the appearance of burial). Mantled 
plains generally have a smooth appearance except for mesas and polygons in 
some areas (figure 4). Polygons seen in figure 4 are approximately the same 
width &s those in bulbous plains (2-6 km) and are outlined by low-albedo 
r.'.rkings or large troughs. 

The gradational relationship of mantled plains to bulbous plains is 
seen in figjre 5. Since craters and other surface features which are high in 
relief are the last features to be obscured with an increase in mantling, 
mantled plains are clearly shown to be overlying bulbous plains. The presence 
of a small bowl-shaped crater in the lower portion of the mosaic In figure 5 
confirms that the obscuring of surface features in the lower portion is due to 



14 

surface mantling and not to cloud cover. 




Dune Deposits 

Dune deposits are characterized by a very low albedo relative to all 
other north polar units and commonly, though not always, by the presence of 
linear and arcuate patterns which have been interpreted as dunes (14). An 
excellent example of the relationship of dunes to mantled and bulbous plains 
Is seen In figure 6. 

Layered Deposits/ Perennial Ice Cap 

Figure 7 shows an example of the layered deposits and the overlying 
perennial ice cap in the north polar region. The reader can refer to numerous 
other papers (8, 13, 14, 15) for detailed description and discussion of these 
units. The present paper will not examine these deposits in detail but 
will instead concentrate on stratigraphic units surrounding the north polar 
cap. 
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Figure 3: Type area for bulbous plains at 72°N, 290°W (Viking frames 
538B01-538B16). Note the characteristic mottled, bulbous appearance 
resulting from bright material within and on the rims of craters and 
surrounding highly-irregular polygons. 
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Figure 5: Viking Orbiter mosaic showing gradation of bulbous plains in the 

W 8EW" "“I? l',‘» heavily-mantled plains In the lower Ji?!; at 
70 N, 260 w (frames 576B10-576B20) . The sharp appearance of the bow I -shape* 
nh^nroi"^^ 6 ? b X arr °^ indicates that the area Is mantled rather than 
obscured by clouds. Also, note the two large “survivor" craters protrudinq 
through bulbous plains In the upper part of the mosaic. 9 






Figure 6: Excellent example of the relationships among bulbous plains, 

mantled plains, and dunes; at 72°N, 59°W (Viking frames 525B01-525B16) . 
There Is a relatively rapid gradation from bulbous plains near the bottom 
of the mosaic (south) to bright mantled plains and dark dunes in the 
north. Dunes In this mosaic are transverse and barchan dunes. Note the 
close, correlation between the presence of dunes and the presence of mantled 
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Figure 7: Viking Orblter frame showing the bright perennial 
layered deposits, which are exposed in the walls of troughs 
(frame 56B61 ) . 
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INTERPRETATION AND DISTRIBUTION OF STRATIGRAPHIC UNITS 

Bulbous and Mantled Plains 
a. Interpretations from visual observations 

Though certainly not the only possible interpretation, bulbous plains 
are interpreted here as having a volcanic origin. This is based primarily on 
its generally low albedo, though a history of near-surface igneous activity 
is also suggested by the presence of what appears to be an exposed dike 
(figure 8). Large "survivor" craters from the underlying older terrain can 
be seen protruding through bulbous plains in figure 8 and figure 5. 

The highly-irregular polygons which characterize bulbous plains are 
somewhat different from polygons seen in lower latitudes. A belt of polygonal 
ground occurs between 40° and 50° N latitude, with typical polygons up to 
20 km across (21). Polygons within bulbous ground are generally not larger 
than 8 km in diameter. Also, the patterned nature of bulbous plains appears 
to be more irregular than the polygonal appearance in the lower latitude 
belt, although this could be a result of the characteristic mottling of bulbous 
plains which tends to accent the irregular polygonal nature. As will be 
shown below, this mottling probably results from bright dust filling the low 
areas surrounding polygons and within craters. Thus, as the amount of 
dust-filling increases, polygon centers can be seen to stand out as highs. 

T he presence of the polygons in the lower latitudes of Mars has been 
attributed to contractive lava cooling or freeze-thaw ice wedging in a 
permafrost layer (21), although a scale problem exists when comparing Martian 
polygons to terrestrial analogs. Polygons formed by ice wedging or lava 
contraction on Earth are typically 100 m wide, whereas those in the 40°-50° 

N latitude belt of Mars are up to 20 km across (21). Although the highly 
irregular polygons within bulbous plains are generally not larger than 8 km. 
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this Is still larger than terrestrial analogs. Still, the Interpretation that 
polygonal ground on Mars may have formed by lava contraction or ice wedging is 
plausible. The unique occurrence of polygonal ground at latitudes greater 
than 40° N tends to favor the interpretation that polygonal ground on Mars 
results from ice wedging (21 ). However, since bulbous plains are confined to 
the northern latitudes, and may be volcanic, it is reasonable that eitier 
process might have formed its characteristic polygons. 

As seen in figure 5, there is a gradation between mantled c * ; ilbous 
plains, such that mantled plains range from lightly-mantled plain* > which 
some features of bulbous plains, such as a relative abundance of small craters 
and a mottled, polygonal appearance, are retained [reference Viking Orbiter 
frames 531B21-531B36], to moderately- and heavily-mantled plains, where 
bulbous characteristics are replaced by a relatively smooth, crater-deficient 
appearance [reference Viking mosaics 211-5573 (frames 71R61-71R62; 71R41-71R45) 
and 211-5562, respectively]. Due to this blanketing nature, mantled plains 
is interpreted as a deposit of wind blown dust covering bulbous plains in 
varying thicknesses. 

Recause of the difficulties involved in defining distinctive boundaries 
within highly gradational units, lightly-mantled plains are henceforth combined 
with bulbous plains to form a distinct mappable unit, (figure 2). lhus. as 
presently defined, bulbous plains include all terrain which retain the 
characteristics of an abundance of small craters and a mottled, bulbous 
appearance. Likewise, moderately- and heavily-mantled plains which as will 
be discussed later could be arbitrarily separated on the basis of relative 
deficiences of small craters, are for convenience combined into a single 
map unit called mantled plains. 

The presence of mesa-like landforms within mantled ground suggests that 
mantled plains have experienced a period of deflation. The origin of the 



Figure 8a: Viking Orbiter mosaic showing what appears to be a linear dike 

and two “survivor" craters within lightly-mantled to moderately-mantled 
bulbous plains; at 70°N, 244°H (frames 505Bol-505B13). Figures 8b and 8c 
are enlargements of these features and allow better discrimination of the 
underlying bulbous nature of the area. 
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Figure 8b: Enlargement of an ancient "survivor" crater seen in figure 8a, 
allowing better discrimination of the underlying bulbous plains character- 
istics, which are somewhat obscured by mantling. 
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polygons within heavily-mantled plains is as poorly understood as the origin of 
the polygons which characterize bulbous plains. The similarity in size, however, 
suggests that either the process of formation is similar, or the polygons of 
mantled plains are controlled by underlying bulbous plains polygons. 

b. Cumulative crater size- frequency distributions 

Cumulative size-frequency distribution plots for small craters ( < 10 km) 
on bulbous and mantled plains are presented in figure 9. The two plots for 
bulbous plains (figures 9a and 9b) represent the areas seen in figures 3 and 
6, respectively. Clusters of similar-sized craters, like those seen in 
figure 3, result from secondary impact of ejecta from large impact events, 
and are generally not representative of the unit as a while. Therefore, all 
craters which appeared to be associated with crater clusters were not included 
in the counts. 

Both plots for bulbous plains display a power-law slope of about -2.0. 

A linear fit to the data gives a one-kilometer intercept "age" (22) of 4000 
craters per 10 km'. 

An unrealistic decrease in the number of small craters counted in 
images can result from imagery degrading as crater diameters approach the 
resolution limits of the images. This results in a shallower slope, or a 
"rolling over," of the cumulative size-frequency distribution curve at small 
diameters. However, Vollover" due to resolution degradation for Viking 
Orbiter imagery occurs at 6 times the pixel resolution (23 ) , or at 1 .7 x 10' 1 km 
for figure 9a. Since in figure 9a the break in slope occurs at 27 times the 
pixel resolution, or at 7.6 x 10~^ km, the slope of -0.7 for smaller crater 
diameters is probably a true feature of the crater population and not a res- 
olution effect. 

Crater size frequency distributions for mantled plains are seen in 



26 


figure 9c, with the plot for heavily-mantled plains representing the area at 
the mouth of Borealis Chasma (figure 4). The plots for both moderately- and 
heavily-mantled plains show a repetition of the shallow -0.7 slope observed 
for smaller craters of bulbous plains. Moderately-mantled plains also display 
a two-segment distribution curve, with a slope of -2.1 for crater diameters 
larger than about 6 km. A slope of -2.1 is equivalent to the -2.0 to -2.2 
slope values for large craters of bulbous plains. 

A composite of figures 9a and 9c is seen in figure 9d. The curve 
segments are for bulbous (curve A), moderately-mantled (curve B), and 
heavily-mantled (curve C) plains. When adjusted for cumulative size-frequency 
functions rather than differential size-frequency functions, the Jones and 
Chapman model for crater obliteration (24) states that two segment curves, 
such as A and B in figure 9d, represent a balance between a production 
function, 

T(x) * f~ C j(x) f T C.(t) dx dt [1] 

* 4 ) a J T-t e 1 

and a cratering-obliteration equilibrium function, 

v(x) e = t c t (T) I °( T ) ^ f a(x) c d U) dx [23 

where x = the crater diameter, t = the retention age of a crater of a 
given diameter, c t (T) = the time dependence of the cratering rate, c^x) = 
the diameter dependence of the cratering rate (c^x) = O' ), o(T) = the 
obliteration rate, and a(x) = the amount of obliteration needed to render a 
crater indiscernable. Since the amount of obliteration, a(x)’ needed to 
remove a small crater is less than that needed for a larger crater, the 
cumulative size-frequency plot for an obliterated (e.g. dust-mantled) surface 
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would show a shallow slope for small diameters, and a steep slope for larger 
diameters. Such is the case for curves in figure 9d. The repetition of 
-2.0 and -0.7 slopes for plots of both bulbous and mantled plains strongly 
suggests that the curves represent an equilibrium, or post- transient, 
cratering-obliteration process. Thus, using the Chapman-Jones model, -2.0 
would define the slope of the production function, while the slope of -0.7 
would represent a cratering-oblitcrat.’on equilibrium function. 

The repetition of the -0.7 equilibrium slope for all crater si 2 o- 

If 

Irequency distribution plots given above implies that the functions D 
and a(x) do not change from one area to another. Otherwise, the slopes 
would not be the same. Thus, making the reasonable assumption that the time- 
dependent cratering rate, c t (T), is also constant for all areas of concern, 
then the value of the equilibrium function, T e (T), for a given diameter 
range, is dependent only on the obliteration rate, o(t). Thus, the ratio 
of equilibrium frequencies, f e , for two such areas, would allow a relative 
obliteration rate, o(t), to be calculated: 


A *(x) ei [ c f (T) 1 / o(T) 2 1 f D b a(x)i c^x)! dx 



For c t (T) != c t (T) 2 , a(x) ! = a(x) 2 , and c d (x)j = c d (x) 2 : 


A y(x) 6i 0(T) 2 

A V(x) ez 0 (T) i 


[ 4 ] 


From figure 9d, the relative obliteration rates for bulbous : moderately 
ma..tled : heavily-mantled plains are calculated as 1:9:88. 
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The change of slope of about 1 between crater-production and obliteration 
equilibrium functions of the north polar units is consistent with the simple 
dust-filling model of Chapman (25). This supports this author's interpretation 
that mantled plains are deposits of wind-blown dust covering bulbous ground 
In varying thicknesses. The break in slope between crater-production and 
obliteration-equilibrium functions is a result of the presence of "survivor 
craters" which have not been obscured by the obliteration process (22). Thus, 
if dust mantling is the only obliterating process acting, the diameter at 
which the function slope changes ( i.e. the "kink point") should be an indi- 
cation of the average thickness of the dust mantle. Based on crater data for 
the Martian surface, Cintala and Mouginis Mark (26) empirically derived that 
the crater depth = 0,104 x diameter®' 9 ^ for craters less than 12.3 km in 
diameter. By substituting the kink-point diameter for moderately-mantled 
plains into the equation, the thickness of the dust layer comprising moderately 
mantled plains is calculated to be about 600 m. Since the crater-function 
kink for heavily-mantled plains does not occur within the range of diameters 
plotted, the thickness of dust deposited in this region is, at least, the 
depth of the largest crater measured, or 770 m. These values are somewhat 
higher than the 200 m maximum predicted by Squyres (15). 

The data points seen on the composite plot (figure 9d) are not data from 
north polar units, but are instead the data from Neukum and Wise (22) for a 
cratered plains unit at Tempe Plateau (30°N, 90°W). The excellent correlation 
between their data and the composite plot for bulbous and moderately-mantled 
plains suggests first that the crater function for bulbous plains may follow 
the moderately-mantled plains curve in the region of diameters greater than 
5 km. Since only smaller craters would be buried by an obliteration episode, 
this is a reasonable assumption if moderately-mantled plains result from 
obliteration of bulbous plains. However, such a theory should certainly be 
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Figure 9a: Cumulative size-frequency distribution plot for small craters 
on bulbous plains. Area represented is from the bulbous plains type area 
seen in figure 3. N = 156 craters; total area = 5,419.1 km . Error bars 
signify one standard deviation confidence intervals and equal ±V7Tfor 
the nth crater. 
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Figure 9c: Cumulative size-frequency distribution plots for small ( $ 13 km) 
craters on mantled plains. The heavily-mantled plains plot represents the 
type area at Borealis Chasma seen In figure 4. N = 9 craters; total area 
« 55,120.6 km 2 . The curve for moderately-mantled plains is for the area at 
76°N, 90°W (Viking frames 70B04, 71B61-71B64). N = 71 craters; total area 
= 24.344.6 km 2 
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Figure 9d: Composite of cumulative crater size-frequency distribution 
curves seen in figures 9a and 9c. Curve A - bulbous plains; curve B - 
moderately-mantled plains; curve C - heavily-mantled plains. Data points 
are from Neukum and Wise (22) for a cratered plains unit at Tempe Plateau, 
30°N, 90°W. 
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tested by plotting the crater distribution curve for greater than 5 km craters 
on bulbous plains. 

Of more Importance Is that the correlation of the bulbous plains 
crater plot with that of the cratered plains unit at Tempe Plateau* places 
bulbous plains within the Martian global stratigraphic sequence as one of the 
relatively old cratered plains units. According to the time scale of 
Neukum and Wise (22), bulbous plains has an "absolute" age of 3.8 b.y. (note: 
according to the same scale, the oldest cratered plains unit, Lunae Planum, 

Is 3.9 b.y.). 

Pune Deposits 

Dunes in the north polar region of Mars are predominantly transverse and 
barchan dunes, although some longitudinal dunes are also present. They are 
similar In size to megadunes (draa) found in the Sahara, Arabia, and the 
Algodones dune fields on Earth (36). 

As suggested below, dune deposits might also occur as "sheet-sand" 
deposits. Figures 1 and 2 show the stratigraphic map and a Mariner 9 mosaic 
of the north polar region of Mars. A visual comparison between the strati- 
graphic map and the Mariner 9 image between 30° - 290° W longitude confirms 
that the dark circumpolar band seen in the mosaic correlates well with the 
presence of dunes. However, in the region 290 o -0°-30° W longitude, the 
dark band extends into much of the area mapped as bulbous plains, whereas 
other areas of bulbous plains, Including the type area, are not as dark. 

Though this region was often cloudy during much of the Viking Orbiter imagery 
sequences, some relatively clear images (e.g. 531B21 - 531B36) confirm the 
duneless, bulbous nature of the area. The presence of the low-albedo band 
In duneless terrain suggests either a very thin "sheet-sand" deposit of dark 
dune material overlying much of the most northern exposures of bulbous plains. 
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or that the albedo of bulbous plains In other areas Is more influenced by 
a thin mantling of bright dust. Notice that even the type area for bulbous 
plains appears to have dust in many depressions (figure 3). 

Actual dunes are generally confined to mantled regions, as illustrated 
in figure 6 and the stratigraphic map in figure 2. If dune deposits overlie 
bulbous plains, they occur predominantly as "sheet sand" deposits rather than 
dunes. This suggests that, unlike bulbous plains, mantled plains provide 
the proper substratum for the accumulation of dunes. Mantling probably 
creates a smooth, crater-deficient surface which might favor dune formation. 

An observed eastward gradation from densely-spaced to more dispersed 
dunes provides valuable information for determining the source of the circum- 
polar dune deposits. Fields in the region of 130°-270° W longitude are exten- 
sive, uninterrupted fields with distinct boundaries and consisting of transverse 
dunes with very narrow interdune areas (figure 10, and Viking Orbiter Mosaics 
211-5269 (rev. 61), 211-5272 (revs. 60,61), 211-5431 (revs. 62,63), and 211-5583 
(rev. 72)). Dunes within the 30°-130° W longitude region occur in more 
dispersed fields and are generally transverse and barchan dunes with wide, 
bare interdune areas (figure 11 and Viking Orbiter Mosaics 211-5269 (rev. 59), 
and 211-5562 (revs. 50,69)). Between 30°-0°-270 o W longitude, dune deposits 
are predominantly in the form of very isolated dunes, or possibly as sheet 
sands. The overall view is one of counterclockwise giadauon frcm dense 
transverse dunes in the region between 270° and 150° W longitude to more 
highly dispersed transverse and barchan dunes and eventually to possible 
sheet sands and isolated dunes. Although the presence of mantling appears to 
be required for dune formation, the degree of mantling does not appear to 
control the morphology of the dune fields. On Earth, transitions from trans- 
verse to barchan dunes and from narrow to wide interdune areas are generally 
related to a decrease in the supply of dune material (27). This relationship 
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Is typically encountered as one moves progressively away from the source of 
sediment. Thus, the situation observed In the north >.olar region may be 
one of a counterclockwise thinning of the supply away from the source. 

Layered Deposits 

The layered deposits are probably a facies of mantled plains In which 
wind-blown dust Is deposited onto the perennial Ice cap rather than directly 
onto the ground ( 15 ). Although It is not clearly understood which changes 
In Mars' orbital parameters most affect the eventual growth and decay of the 
perennial Ice caps (3), the extent of layered deposits at the poles probably 
Indicates the maximum extent of the perennial ice caps in the past. 

Relationships of units at Borealis Chasma may Illustrate an inter fingering 
of mantled plains and layered deposits. Layering can be seen In the scarp 
which overlies mantled plains at the "mouth" of Borealis Chasma (figure 4). 
However, the surface of much of Borealis Chasma does not display the laminations 
within minor topography as Is typical of north polar layered deposits. The 
surface of Borealis Chasma instead has the appearance of mantled plains, 
suggesting a transition from mantled plains to layered deposits and back to 
mantled plains. Such a transition is not unreasonable since the extent 
of the perennial Ice cap is expected to fluctuate with time. 

Comparison with Previous Research 

The north polar stratigraphy presented in this paper has been compared 
to the stratigraphy proposed by Soderblom et. al. (11), and that of Squyres 
(15). Although the preliminary stratigraphic map of Soderblom et. al. has 
been greatly altered, their basic stratigraphic sequence for the north polar 
region has not changed drastically. Bulbous plains are certainly a sub-unit 
of the "mottled cratered plains" unit which includes a large portion of the 
northern hemisphere. However, bulbous plains with their characteristic highly- 
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Irregular polygonal appearance are probably confined to the northern portion 
of the mottled cratered plains. Mantled plains are equivalent to the debris 
mantle and etch-pitted plains In the stratigraphic sequence of Soderblom et. 
al., whereas dune deposits replace the rippled plains. 

Although the stratigraphic sequence proposed In this paper agrees In 
theory with that proposed by Squyres (15» 28). there are very serious discrep- 
ancies In the defining and mapping of the cratered plains and mantled units. 

The previous author's type locale and only exposure of his underlying "cratered 
plains" unit Is at the mouth of Borealis Chasma. This same area, seer. In 
figure 4 and represented by the crater plot for heavily-mantled plains In 
figures 9c and 9d, is the present author’s most heavily-mantled unit In the 
north polar region. Both the imagery and the crater size-frequency distribution 
data for this area show an obvious deficiency of craters relative to moderately- 
mantled and bulbous plains. Also, as discussed earlier, the shallow slope of 
-0.7 confirms that the unit is a highly-obliterated surface, probably resulting 
from a dust mantling which may be at least about 700 m thick. Units repre- 
sented by curves A and 3 in figure 9d are one to two orders of magnitude 
more cratered than the "cratered plains" unit of Squyres. Both units repre- 
sented by curves A and B have been mapped by the previous author as a debris 
mantle which formed by dust obliteration of craters on the "cratered plains" 
unit represented by curve C. Such an Interpretation of these units Is 
obviously not consistent with the crater size-frequency distribution function 
for these areas. 
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Numerous investigators (29, 30, 31 ), including the present author 
(this paper), have independently inferred near-surface wind directions 
within the north polar region of Mars, on the basis of dune orientations. 
Figure 12 shows dune orientations mapped from orthographic images which 
minimize distortions resulting from "oblique viewing" of the spacecraft. 

The map presents data in the form of orientations of arcuate (barchan) 
dunes and the trend of linear (transverse and longitudinal) dunes. The 
crisscross pattern cf the dunes between 150° and 210° W longitude can be 
seen in figure 10. 

Caution must be exercised when interpreting wind directions based on 
dune orientations derived from Viking Orbiter images. In such cases, barchan 
forms provide the most reliable indicators of near-surface wind orientations. 
Transverse and longitudinal dunes in the circumpolar images generally only 
allow determination of a wind trend rather than a wind direction. Thus, 
interpretations based on transverse and longitudinal dunes must rely on 
indirect evidence, such as wind streaks or directions based on nearby barchan 
dunes. Tsoar et. al. (31) have attempted to determine absolute wind dir- 
directions based on the asymmetry of transverse dunes. However, it is again 
cautioned that the resolution of Viking north polar images may not allow 
reliable separation of true dune slope asymmetry from illusory effects dun 
to shadows or oblique viewing. 

In figure 13 is the present author's interpretation of wind directions 
based primarily on orientations shown in figure 12. Interpretations of wind 
directions for the transverse dune fields between 120° and 240° W longitude 
relied on the assumption of a continuity of wind flow relative to barchan 
dunes which occur at the edge of the field. Wind directions in figure 13 
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indicate a predominantly counterclockwise flow of dune-driving winds around the 
ice cap along with dune modifying winds which spiral clockwise off of the ice 
cap, possibly during a change in season. Such an interpretation of wind 
directions is consistent with north polar wind circulation patterns derived 
from streak orientations (29), from lee cloud orientations (32; R.G. French, 
personal communication, 9/25/79), and from other interpretations based on 
dune orientations (29 , 30, 31 ). Leach (30) attributes the clockwise pattern 
to the influence of Coriolis force acting on katabatic winds spiralling off 
the ice cap. 

Thus, data derived from dune, wind streak, and lee cloud orientations 
indicate a polar atmospheric circulation pattern consisting of winds travelling 
around the ice cap in a counterclockwise direction, modified by katabatic 
winds spiralling clockwise off of the ice Cup. Such a pattern is similar to 
the wind circulation around the Antarctic ice cap on Earth (30). 
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Figure 12: Dune trends within the north polar region of Mars. (1) linear 
trends within dune fields, and grooves in ice, (2) direction of arcuate forms 
in dunes (asymmetry of symbols corresponds to asynmetry observed), (3) complex 
dunes with no distinct linear or arcuate patterns, and (4) possible sheet sand 
deposits. Symbols signify general trends of dunes rather than individual dunes. 
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IMPLICATIONS FOR THE SOURCE OF CIRCUMPOLAR DUNE MATERIAL 

As discussed earlier, the gradation of extensive, uninterrupted 
transverse dune fields to scant fields with widely spaced transverse and 
barchan dunes implies a general decrease in the supply of dune material as 
one moves eastwardly (or counterclockwise) away from the area between 270° 
and 180° W longitude. Circumpolar winds were found to flow predominantly 
counterclockwise, suggesting a general thinning of dune deposits in a downwind 
direction. This, combined with the presence of extensive exposures of 
bulbous plains upwind of extensive dune fields, may indicate that dune material 
is stripped from bulbous plains and deposited on mantled regions. 

Such a conclusion might be further supported if one considers the dark 
circumpolar band extended into duneless bulbous plains areas. Two possible 
interpretations were suggested earlier to account for this, one being that 
a very thin sheet sand deposit of dark dune material blankets bulbous plains 
in this area, and the other suggesting that bulbous plains in this area are 
more completely stripped of a thin bright dust layer than bulbous plains in 
other areas. The first interpretation is consistent with dune material being 
continually derived from bulbous plains only to be slowly stripped from the 
region and deposited onto mantled plains downwind. If the second interpretation 
were true it would lend very strong support to the suggestion that bulbous plains 
are the source of the dune deposits. Not only would it provide evidence that 
this possible source area is significantly stripped, but more importantly 
it would indicate that bulbous plains and dune deposits have the same low 
albedo. 
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CONCLUSION 

The stratigraphic sequence of the north polar region of Mars can be 
separated into four distinct informal units: (a) bulbous plains, (b) mantled 
plains, (c) dune deposits, and (d) layered deposits/perennial Ice. 

Bulbous plains is a relatively old "mottled" cratered plains unit 
equivalent in crater-distribution age to the cratered plains unit at Tempe 
Plateau (90°W, 30°N). According to the time scale of Neukum and Wise, the 
"absolute age" of bulbous plains is 3.8 b.y.. Due to its relatively low 
albedo and because of near-surface volcanics suggested by the presence of 
what appears to be a dike, bulbous plains are interpreted as having a volcanic 
origin. Polygons within bulbous plains are 3-8 km and may have formed by 
contractive cooling of lava or freeze-thaw ice wedging within a permafrost 
layer. 

Mantled plains are interpreted as wind-blown dust deposits covering 
bulbous plains in varying thicknesses. A deficiency of small craters seen 
in crater size-frequency distribution functions for bulbous and mantled 
plains supports the interpretation that mantled plains formed by obliteration 
of bulbous plains, while a change of slope of about 1 between crater-production 
and obliteration-equilibrium functions for north polar units is consistent 
with obliteration by dust-filling. The relative surface-obliteration rates 
are 1:9:88 for bulbous, moderately-mantled, and heavily-mantled plains, 
respectively. 

Previous stratigraphic mapping proposed by Squyres for the north polar 
region of Mars is inconsistent with interpretations based on crater size- 
frequency distribution data. Such evidence strongly supports that his 
"cratered plains" unit at the mouth of Borealis Chasma is instead a heavily- 
mantled plains, whereas many areas mapped by the previous author as "debris 
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mantle" are instead cratered plains. 

Layered deposits probably represent a facies of mantled plains in which 
dust is deposited onto the perennial ice cap rather than directly onto the 
ground. Thus, the extent of the layered deposits at the north and south 
poles may indicate the maximum extent of the perennial ice caps in the past. 

A possible interfingering of layered deposits and mantled plains can be 
seen at the mouth of Borealis Chasma. 

Dunes in the north polar region of Mars have a very low albedo and 
occur as transverse, barchan, and to a minor extent, as longitudinal dunes. 
Actual dunes are generally confined to mantled plains, suggesting that 
mantling provides a proper substratum for dune formation. The extension of 
the dark circumpolar band into duneless areas mapped as bulbous plains 
suggests that thin sheet sand deposits overlie bulbous plains in this area 
or that the most northern exposures of bulbous plains have been more effect- 
ively stripped of bright dust than bulbous plains in other areas. If the 
latter interpretation is true, it indicates that the albedo of bulbous 
plains is similar to the albedo of the dunes. The distribution of dune forms 
and dune field styles reveals thinning of the supply of dune material in an 
eastward direction away from the area between 270° and 180° W longitude. 

As determined from dune orientations, near-surface north polar winds 
consist of a counterclockwise flow around the ice cap modified by katabatic 
winds spiralling clockwise off the cap. This is analogous to the atmospheric 
circulation pattern around the Antarctic ice cap on Earth. 

The results above suggest a general thinning of the supply of dune 
material downwind of the area between 270° and 180° W longitude. Immediately 
upwind of this area are extensive exposures of bulbous plains. Thus, the author 
suggests that dune material is being stripped from bulbous plains and is 
accumulating as dunes downwind on mantled plains. 
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RiTROrUCiMUTY OF THE 
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mm*tM***m*m******tttt***m*M**m***t«*t*******t****«*i 

SUBROUTINE CRhODl *FTK (CRATER DISTRIBUTION NODULE «1> 

MM************************************************************* 


»»»> CUMULATIVE SIZE-FREQUENCY DISTRIBUTION PLOT «««< 


HIRE BOTTS 2/U/79 


COMPLETE SUBROUTINE FOR CALCULATING AND PLOTTING CUMULATIVE 
8 I 2E- FREQUENCY DISTRIBUTION DATA (ALONG WITH ERROR BARS)* 

OUTPUT INCLUDES TABLE OF VALUES NEEDED FOR PLOTTING CUMULATIVE 
SIZE-FREC. DISTRIBUTION# AND GRAPH OF VALUES AND ERROR LIMITS. 

INPUT I D « ARRAY OF UNBINNED CRATER DIAMETERS ARRANGED IN 
DESCENDING ORDER 
N - NUMBER OF CRATERS 

NDIM * ROW DIMENSION GIVEN 'CUMPLT' IN MAIN PROGRAM 
(DIMENSION GIVEN 'BARCUM' SHOULD BE 2*NDIM> 

A ■ TOTAL AREA <KM**2) OF CRATER COUNTS 
XAXXS1 • LOWER BOUND OF X-AXIS 
XAXIS2 » UPPER BOUND Or X-AXIS 
YAXIS1 « LOWER BOUND OF Y-AXIS 
YAXX82 * UPPER BOUND OF Y-AXIS 

NOTES I IF LIMITS OF GRAPH ARE NOT KNOWN# SET XAXIS1 AND XAXIS2 
EQUAL TO 0.0 .PROGRAM PLOT • FTN WILL THEN COMPUTE 
ITS OWN INCREMENTS. 

X INCREMENT « (XAXIS2 - XAXISD/100. 

Y INCREMENT - < YAXIS2 - YAXISD/50 . 

(PROGRAM WILL ASK IF YOU DESIRE A REPLOT) • 

RETURNS $ 

DXA - DIAMETERS OF CRATERS 

CUM « CUMULATIVE NUMBER OF CRATERS LARGER OR EQUAL TO 
DIAMETER AT THAT POINT 
DIALOG « LOG (D) 

CUMLOG - LOG (CUM) 

ERRLOG « LOG (1ST SET OF ERROR LIMITS) 

AIRLOG - LOG (2ND SET OF ERROR LIMITS) 

CUMPLT - 2-COLUMN ARRAY OF MEASURED DATA WHERE I 
CUMPLT ( I # 1 ) » DIALOG < I ) 

CUMPLT < I #2) » CUMLOG ( I ) 

BARCUM « 2-COLUMN ARRAY OF ERROR BAR DATA WHERE: 

" BARCUM < I # 1 ) » DIALOG ( I ) 

BARCUM ( I # 2 ) » ERRLOG ( I > 

BARCUM ( I+i » 1 ) * DIALOG ( I ) 

BARCUM ( 1+1 #2) * AIRLOG(I) 

J - DIMENSION VARIABLE OF DIALOG# CUMLOG# ERRLOG# AIRLOG 


SUBROUTINES NEEDED: 

CUMFRQ.FTN (CALCULATES VALUES NEEDED FOR PLOT) 


SQUEEZ.FTN (PREPARES MEASURED DATA FOR PLOTTING PROGRAMS) 
BOUISH.FTN (PREPARES ERROR BAR DATA FOR PLOTTING PROGRAMS 
PLOT. FTN (PLOTS DATA + ERROR BARS) 


SISTER PROGRAMS! 

CRM0D2.FTN (DIFFERENTIAL SIZE-FREQUENCY PLOT) 

CRM0D3.FTN (BINNED RELATIVE SI. E-FREQUENCY PLOT) 


f 



c 


SUBROUTINE CRNOD1 (DtNrNQIMr AtXAXISl .XAXIS2. YAXIS1 r YAXIS2?DIA»CUH ( r 
1 DIALOG r CUMLOG r FRRI 00 . A I RLOG r CUMPL T . BARCUM r J > 

DIMENSION D(N) *DIA(N) • DIALOG (N) r CUM < N ) * CUMLOG ( N ) r ERRLOG(N) r AIRLOG(N ) 9 
tCUMPLKNDIM* 2> rBARCUM(24NDXMr 2) 

C * 

C CALL ROUTINE FOR CALCULATING CUNS 

C 

CALL CUNFRO (DrNrAtDIArDI ALOB 9 CUM f CUMLOO 9 ERRLOG r AIRLOG r J ) 

C 

C LIST VALUES NEEDED FOR PLOTS D 9 CUM » D I ALOG 9 CUMLOG r ERRLOT » AIRLOG 

C 

WRITE (Sr 100) 

WRITE (Sr 120) 

URITE (Sr 140) 

WRITE (5r ISO) 

100 FORMAT ('l'r 72 ('*')) 

120 FORMAT ( '0' rSXr 'CUMUL ATIVE SIZE-FREQUENCY DISTRIBUTION' > 

140 FORMAT ( '0' r72( '* ' ) > 

ISO FORMAT ('0 DIAMETER' rSXr 'GUM. NO. /AREA' tSXr 'LOG (DIAMETER) ' 

1 r 5Xt 'LOG (CUM. NO. /AREA) ' rSXr 'LOO (ERROR LIMITS)') 

DO 200 K-l r J 

WRITE (5r 180) DIA(K) »CUM<K) .DIALOG(K) rCUMLOG(K) rERRLOG(K) r AIRLOO(K) 
180 FORMAT (' ' rF10. 3.8X.G10. 3r4X. ' ! ' .4X.F10.3r 13Xr F10. 3r 3X»2F10.3) 

IF (AIRLOG(K).EQ. -9999*999) AIRLOG(K) - CUMLOG(K) 

200 CONTINUE 

C 

C PREPARE DATA FOR PLOT 

C 

CALL SPUEEZ (DIALOG r CUMLOG t J r NDIM 9 CUMPLT ) 

CALL SQUISH ( DIALOG r ERRLOG r AIRLOG 9 J r 24NDIM r BARCUM r J J > 

C 

C PLOT DATA AND ERROR LIMITS 

C 

URITE (Sr 300) 

WRITE (Sr320) 

WRITE (5r 330) 

300 FORMAT ('OX • LOG( DIAMETER (KM) >') 

320 FORMAT (' Y « LOG( CUMULATIVE NUMBER OF CRATERS PER KM**2 )') 

330 FORMAT (' PHI » ( CUM. NO. +- SQRT( CUM. NO. > ) / AREA') 

CALL PLOT (CUMPLT r JrNDIMr BARCUMr JJ*24NDIMr XAXI81 rXAXIS2 r YAXIS1 r YAXIS2) 

C . ... 

C ASK FOR REPLOT OF GRAPH 

C 

490 WRITE (SrSOO) 

SOO FORMAT ( 'tUOULD YOU LIKE TO REPLOT THIS GRAPH? ( O-NOr 1-YES) S ' ) 

READ (5r50S) NSER 
SOS FORMAT (ID 

IF (NSER.EQ.O) GO TO 600 
URITE (Sr 520; 

323 FORMAT ('tLIST XMIN. XMAX. YMINr YMAX 5 ') 

READ (5.525) XAX1S1 .XAXIS2.YAXIS1. V AXIS2 
323 FORMAT (4F10.5) 

CALL PLOT(CUMPLTr Jr NDIM. BARCIJMr JJr 2*NDIMr XAXIS1 r XAXIS2. YAXISt . YAXIS2) 
GO TO 490 
600 RETURN 

END 
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SUBROUTINE CUMFRO.FTN (CALCULATES CUM- SIZE FREQUENCY VALUES) 


NIKE S0TT8 1/31/79 

CALCULATES VALUES NEEDED TO PLOT CUMULATIVE SIZE-FREQUENCY 
DISTRIBUTIONS ACCORDING TO STANDARDS SUGGESTED BY 
•STANDARD TECHNIQUES FOR PRESENTATION AND ANALYSIS OF 
CRATER SIZE-FREQUENCY DATA* (NASA TM-79730). 

INPUTS 

D(N) * ARRAY OF MEASURED DIAMETERS (KM) OF N CRATERS 
A - AREA (KM**2> OF TOTAL PICTURES COUNTED 

OUTPUT! 

DIACJ) • ARRAY OF DIAMETERS HAVING DIFFERENT VALUES 
DXALOO(J) - LOG< DIA<J> ) 

CUM(J) » CUMULATIVE NO. OF CRATERS LARGER OR EQUAL TQ 
DIA(J) • 

CRRLOO(J) AND AIRLOG(J) * LOG ( PHZ(J) ) 

PHK J) * CCUM( J) ♦- SQRT< CUM(J) )3 / A 
J • TOTAL NO. OF CRATERS WITH DIFFERENT DIAMETERS 

note: array d must list diameters in decendxng order from 
0(1) TO D(N). use subroutine sortc. 

NEGATIVE INFINITY MARKER « -9999.999 . 


SUBROUTINE CUMFRQ<D»Nf A tDIAr DIALOG* CUM tCUMLQGf ERRLOG* AIRLOG, J) 
DIMENSION D(N) * DIA<N) t DIALOG(N) rCUM(N) rCUMLOG(N) vERRLOG(N) »AIRL0G(N) 

M«0 

DO 100 1*1 »N 
IF (I.EQ.N) 60 TO 60 
IF (P(I> »EQ.D(I+1) ) 00 TO 100 
M«M4i 

CALCULATE VALUES FOR EACH DIFFERENT DIAMETER 
DIA(M) - D(I> 

DIALOO(M) » AL0G10(D(I) > 

C - 1.0 * FLOAT < I ) 

CUN(M) - C/A 

CUMLOO(M) - AL0610(CUM(M> > 

CRRLOQ(M) - AL0G10( (C ♦ SQRT(C)>/A> 

IF (I.OT.I) GO TO 90 

AIRLOG(M) - -9999.999 

GO TO 100 

90 AIRLOO(M) - ALOGlOdC- SORT(C)>/A) 

100 CONTINUE 

J»H 

RETURN 

END 
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c ******************************************************************* 

c 

C SUBROUTINE SORTC*FTN (SORTS COLUMNS OF AN ARRAY) 

C 

C ******************************************************************* 

c 

C MIKE POTTS 1/31/7? 

C 

C SORTS ELEMENTS OF EACH COLUMN OF ARRAY A(NrM) INTO DESCENDING ORDER 

C FROM TOP TO BOTTOM AND RETURNS AS APRIM£<N>M) 

C 

C N • ACTUAL NUMBER OF ROUS IN ARRAY A 

C M - ACTUAL NUMBER OF COLUMNS IN ARRAY A 

C NPXH - ROW DIMENSION GIVEN ARRAY A IN MAIN PROGRAM 

C 

c 

SUBROUTINE SORTC ( A » N t ND I M » M r APR I ME ) 

DIMENSION A(NDIMtM) f APRIME(NDIMfN) 

DO 20 I«1>N 
DO 20 J*l t M 

APRINE ( I » J) - A(X*J) * 

20 CONTINUE 

L-N-l 

DO 100 K*1 t L 
DO 100 J»1*M 
DO 100 1*1 t N 
ZF(I*EQ«N) 00 TO 100 

IF<APRlME(It J).LE.APRIME(I+1*J>> GO TO 50 
GO TO 100 

50 TEMP - APRXMECXr J) 

APRINE ( Xf J) * APRIME< 1+1 » J> 

APRXME< 1+1 t J) - TEMP 
100 CONTINUE 

RETURN 
END 


C ***************************************************************** 

c 

C SUBROUTINE SQUEEZ.FTN (SQUEEZES 2 ARRAYS INTO 1 ARRAY) 

C 

c **************************************************************** 

c 

C MIKE BOTTS 2/1/79 

C 

C TURNS 2 ONE-COLUMN ARRAYS INTO 1 TWO-COLUMN ARRAY. 

C XY( I • 1 ) « XU> 

C XY(I>2) - Y(I) 

C 

C I ■ ACTUAL NUMBER OF X AND Y VALUES 

C NDIM - ROW DIMENSION GIVEN XY IN MAIN PROGRAM 

C 

C GOOD FOR PREPARING DATA FOR SUBROUTINE PLOTFR.FTN* 

C 

c 

SUBROUTINE 5GUEEZ(X> Y» I rNDIMrXY) 

DIMENSION X(I)f Y ( I ) » XY<NflIM>2) 

DO 100 J*1'I 
XY(Jrl) - X<J> 

XY(Jt2) ■ Y( J) 

100 CONTINUE 

RETURN 
END 
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SUBROUTINE SQUISH. FTN (SQUISHES 2 UNEQUAL ARRAYS INTO X ARRAY) 
**t*»tt«M********M*«****«*«*«*M*«**»*******«t******t***««*tt*« 


HIKE S0TT8 2/2/79 

SPECIAL ROUTINE NEEDED TO PUT ERROR BAR DATA INTO A 2-COLUMN 
REQUIRED BY SUBROUTINE PLOTER.FTN. 


VALUE! J) ■ VALUE PLOTTED ON X AXIS (SANE AS X VALUE FOR REAL 
DATA POINTS). 

ERRKJ) - ONE SET OF ERROR BAR LIMITS 

ERR2!J> - OTHER SET OF ERROR BAR LIMITS 

J - ACTUAL NUMBER OF 'VALUE' POINTS 

NDIM - ROW DIMENSION GIVEN 'BAR' IN MAIN PROGRAM 


RETURNS BAR< JJ»2) WHERE JJ • J*2 AND! 
BAR< 1*1) - VALUE < 1 ) « 

BAR<2tl> - VALUE ( 1 > 

BAR(3tX> - VALUE (2> 

BAR<4f 1 ) » VALUE! 2 > 


BARd »2) * ERR1!1> 
BAR(2»2) « ERR2 ( 1 > 
BAR!3»2) * ERR1<2) 
BAR(4»2) - ERR2<2) 


BAR! JJtl ) - VALUE! J) 


BAR! JJ»2> - ERR2! J) . 


USED BY MAIN PROGRAM CRPLOT.FTN (CRATER PLOT) WHICH PRODUCES 
CRATER SIZE-FREQUENCY PLOTS WITH ERROR BARS. 


SUBROUTINE SQUISH!VALUEp£RR1pERR2p J tNDIMrBARt JJ) 
DIMENSION VALUE! J) »ERR1< J) pERR2( J) p BAR < NDIM f2> 

JJ » Jt2 
N-0 

DO XOO I«1pJ 
N*N41 

BAR(NpI) - VALUE ! I ) 

BAR!Np2) » ERRKI) 

N*Nf 1 

BAR(NpI) - VALUE! I) 

BAR!Np2> » ERR?! I ) 

ISO CONTINUE 
RETURN 
END 
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(PLOTS 2 SETS OF SATA ON ONE GRAPH) 


SUBROUTINE PLOT.FTN 


HIKE BOTTS 2/16/79 

MODIFIED VERSION OF DAVIS' PLOTER.FTN «: GRAPH BOUNDARIES CAN NOW 
BE GIVEN BY USER OR CALCULATED BY PROGRAM. IF CALCULATED BY 
PROGRAM# ROUTINE NOT GUARANTEED TO GIVE ‘CLEAN* INTERVALS. 

• (DAVIS# J.C. (1973)# 'STATISTICS S DATA ANALYSIS IN GEOLOGY' #PG 211 

PLOTS TWO VARIABLES ON ONE GRAPH. IF THE TWO VARIABLES CONSIST OF 
DATA + ERROR LIMITS# IT IS SUGGESTED THAT THE DATA ARRAY » Ql# 

AND THE ERROR ARRAY - 02. 

INPUTS 

Ql • FIRST SET OF VARIABLES TO BE PLOTTED# WHERE 
01(1# 1 > • VALUE PLOTTED IN X DIRECTION 
01 (I#2) - VALUE P. 3TTED IN Y DIRECTION 
02 - SECOND SET OF VARIABLES TO BE PLOTTED 
N1 » ACTUAL NUMBER OF ROWS OF 01 

N2 • ACTUAL NUMBER OF ROUS OF 02 (SET » 0 IF ONLY ONE DATA 
SET IS TO BE PLOTTED) 

N3 * ROW DIMENSION QIVEN 01 XN MAIN PROGRAM 

N4 - ROW DIMENSION GIVEN 02 IN MAIN PROGRAM 

XAXI81 • LOWER BOUND OF X AXIS 

XAXXS2 - UPPER BOUND OF X AXIS 

YAXIS1 - LOWER BOUND OF Y AXIS 

YAXIE2 * UPPER BOUND OF Y AXIS 

NOTES! IF YOU WISH BOUNDARIES CALCULATED FOR YOU f SET XAXIS1 AND 
XAXIS2 EQUAL TO 0.0 . 

IF YOU WISH TO GIVE YOUR OWN BOUNDARIES# NOTE THAT : 

X INCREMENT « (XAXIS2 - XAXIS1 1/100 
Y INCREMENT - < YAXIS2 - YAXISD/50 . 


OUTPUT! 

* - 18T DATA SET 
- * 2ND DATA SET 

0 * OVERLAPPING POINT OF 1ST S 2ND SETS 


SUBROUTINE PL0T(Q1 »N1 #N3t 02tN2#N4# XAXIS1 # XAXIS2# YAXIS1 #YAXIS2> 
DIMENSION Q1(N3»2) #02<N4»2) # XOUT ( 101 ) #XX( 11 ) 

DATA XBLNK/' ' /fll/'I V# IPLUS/' + V# IBAR/'_'/» IMINUS/'-V# TSTAR/'#'/ 
DATA IO/'O'/ 

IF XAX2S1 AND XAXIS2 • 0.0 CALCULATE BOUNDARIES AND INCREMENTS 

NX1 » INT(XAXISl) 

NX2 ■ INT (XAXIS2) 

IF ( NX 1 » NE * 0 • OR • NX2 « N£ » 0 ) GO TO 100 


IF BOUNDARIES NOT GIVEN# DETERMINE MAX « MIN 

FIND MAXS AND MINS OF DATA SET 01 

XAXXS1 « Ql(l#l> 

XAXIS2 - XAXI81 
YAXZS1 - Ql < 1#2) 

YAXI82 ■ YAXIS1 
DO 90 J-1#N1 

IF (Q1(J#1) .LT.XAXX81) XAXIS1»Q1 < J# 1 > 


IF <01 < JfI ) *GT*XAXIS 2) XAXIS2-01 < Jtl > 

XF <Q1 < J»2) «LT«YAXIS1 ) YAXIS1-01 < Jf2> 

XF (Q1 ( Jr2) *GT* YAXI32) YAXIS2-01 < Jf2> 

•0 CONTINUE 

XF <N2*EQ«0) 00 TO 100 
C 

C FIND MAX8 AND NXNS OF DATA SET 02 

C 

DO 70 J-t»N2 

XF <Q2( J» 1 ) *LT*XAXIS1 ) XAXIS1«Q2< J»l> 

XF (02<Jf1)»0T *XAXIS2) XAXIS2-Q2< JfI > 

XF (Q2( J»2) *LT* YAXIS1) YAXIS1*02< Jf2> 

IF (Q2< J»2) *UT* YAXIS2) YAXI8?»Q2< Jf2> 

70 CONTINUE 

C * 

C CALCULATE INCREMENTS 

C 

100 DX - (XAXIS2 - XAXI6D/100* 

DY • (YAXI82 - YAXI81 >/50* 

C 

C OCT ORAPH LXNC8 AND DIVISIONS SET UP 

C 

WRITE <5f103> 

10S FORMAT < '0Y-AXX8' »//) 

Y • YAXI82 
DO 500 I-ltSl 

XF <H0D<I-I»5>.EQ.0> 80 TO 120 
DO 110 J«1f101 
XOUT(J) - I BUNK 

XF (M0D<J~1f10) *EQ.O) IOUT(J)«XI 
110 CONTINUE 

00 TO 140 

120 DO 130 J»1 v 101 

XOUT(J) ■ I MINUS 

XF 4H0D( J**lt 10) »£Q«0) IOUT(J) - IPLUS 
130 CONTINUE 

C * 

C PREPARE DATA SET 01 FOR PLOTTING 

C 

140 DO 150 J*1»N1 

IY - NINT( <Q1( J«2> - YAXIS1 )/ DY ) 4 1 
XF i <52-1 > *NE. I Y) GO TO 150 
XX ■ NXNT( <Q1( JfI ) - XAXIS1 )/ DX ) 4 l‘ 
XOUT(XX) - I8TAR 
150 CONTINUE 

C 

0 PREPARE DATA SET 02 FOR PLOTTING 

c 

XF <N2.L£*0> GO TO 200 
DO 140 J*1 tN2 

XY - NXNTC (Q2C J*2> - YAXIS1)/ DY ) 4 1 
XF < <52-1 > *NE* IY) GO TO 140 
IX - NINT ( (Q2< JfI ) - XAXIS1 )/ DX ) 4 1 
XF ( IX*GT • 100) IX » 100 
IF (lOUT(IX).EO.IO) 00 TO 140 
IF <XOUT < IX) *EQ«ISTAR) GO TO 155 
I OUT (IX) - IBAR 
GO TO 140 

155 XOUT < IX) - 10 

140 CONTINUE 

C Ml 

C PLOT ONE LINE OF PLOT 

C 

200 XF <M0DU**1»5).NE'0> 80 TO 300 

WRITE <5f 250) YfIOUT 
250 FORMAT <1X»F10.5»1Xf101A1> 
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Y - Y-5.06DY 
00 TO 500 

300 WRITE (5*350) IOUT 

350 FORMAT <12X*101A1> 

500 CONTINUE 

C 

C • PRINT BOTTOM LEGEND 
C 

XXX * XAXIS1 
DO 600 1*1 r 1 1 
XX(I> * XXX 
XXX - XXX + 10.06DX 
600 CONTINUE 

WRITE (5*700) (XX<I)*W*11.2> 
700 FORMAT <8X*6<F10 . 5 * 10X) ) 

WRITE <5*000) (XX( I ) * 1=2* 10*2) 
800 FORMAT < 10X*5(F1O . 3 * 10X) ) 

WRITE (5*900) 

?00 FORMAT ( SOX * 'X-AXIS') 

RETURN 

END 



appendix c 

Crater Size-Frequency Data 

u: prxeuam' ity of the 

„ . . ORIGINAL PAGE K POOR 

Bulbous Ground i 72 N # 290* w 

framea uaod » 538B03, 5381*05, 538308, 538B09, 538BI2 
total area , 5*19.10? km 2 
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number of craters * 1 56 

Alamcter (km) 


Cumulative no.^. (km” 2 ) 


19.968 

.133 * 10"3 

5.233 

.266 

2.319 

.399 

2.093 

.532 

2.092 

.665 

2.000 

.798 

1.9*3 

.931 

1.5*2 

.106 X 10” 2 

l.*95 

• 133 

l.*69 

.1*6 

1.395 

.160 

1.3*5 

.186 

1.306 

.200 

1.2*8 

• 213 

1.196 

.2*0 

1.175 

. 266 

1.160 

.279 

1.121 

.319 

1.101 

.333 

1.077 

• 373 

1.0*6 

.399 

1.015 

.*39 

1.000 

.*52 

.998 

.*66 

.972 

.5*6 

.955 

.599 

ON 

CO 

.652 

.681 

.719 

.870 

.732 

.858 

.7*5 

.822 

.785 

.797 

.798 

.780 

.812 

.768 

.825 

• 

09 

.905 

.73* 

.971 

.725 

.102 X 10' 1 

.702 

.106 

.692 

.108 

.691 

.113 

.673 

.11* 


' ? 
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(continuod) 


diameter (km) 

cumulative no./ari 

• 661 

.121 X 10". 1 

.652 

• 124 

.624 

.125 

• 615 

.130 

,614 

.132 

.609 

.133 

.598 

.142 

• 587 

.150 

.580 

.157 

.538 

.162 

.523 

.164 

• 51* 

.166 

.507 

.169 

.468 

.170 

.462 

.17? 

.461 

.176 

.449 

,188 

•441 

.196 

.*35 

.200 

.390 

.204 

.385 

.206 

.384 

.209 

.374 

.213 

.367 

.214 

.362 

.216 

.312 

.228 

.308 

.229 

.307 

.237 

.299 

.244 

.294 

.250 

.290 

.253 

.234 

,262 

.231 

.267 

.230 

.273 

.224 

.279 

.220 

.283 

.217 

.285 

.159 

.286 

.156 

.294 

.154 

.315 

.149 

.329 

.147 

.33* 

.1*5 

.337 

.077 

.341 

.075 

.3*3 

.072 

.3*5 


(km“ 2 ) 
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Mtow flrsana • 7**n, 59*w , . 

, i ■ * ' • • ’ i ■ t • - 

frtoos upod i 525801 - 525BO6 1" M.’*' n^jR 

total «r«* 1 5400.687 Km* * ‘ ’ ' 

number of orators 1 97 


diameter (km) 

cumulative 

3.166 

.185 : 

2.174 

.555 

1.871 

.741 

1.594 

.130 : 

1.295 

.146 

1.160 

«A 

CO 

• 

1.151 

. 222 

1.015 

.259 

1.00? 

.296 

,870 

.352 

.725 

.444 

.720 

.537 

.580 

.574 

.576 

.778 

.435 

.107 1 

CM 

• 

.122 

• 290 

.137 

.288 

.156 

.145 

.167 

.144 

.180 



62 


Modorptoly-wftnfl^JLrQund 

1 ?6 *n, 9 <J*w 

frames used i 70B04, 

71B61 - 71B64 

total area i 24,344.504 km 2 

number of craters i 

n 

diameter (km) 

cumulative no. /area 

12.878 

.411 X 10” 1 * 

10.231 

.822 

9.625 

.123 X 10" 3 

6.648 

.164 

6.576 

.205 

6,241 

.246 

5.540 

.298 

4.818 

• 329 

3.942 

.370 

3.042 

.411 

2.847 

'.452 

2.750 

.493 

2.466 

.534 

2.200 

.575 

2.063 

.616 

1.787 

.657 

1.650 

.698 

1.314 

.822 

1.095 

.904 

.962 

.986 

.876 

.107 X 10“ 2 

.831 

.111 

• 

00 

u 

0 

.115 

.825 

.119 

.766 

.127 

.687 

.131 

.657 

.152 

.550 

.156 

.548 

.168 

.438, 

.177 

.416 

.181 

.411 

.185 

.329 

.189 

.277 

.205 

.275 

.209 

.274 

.230 

.219 

.271 

.138 

.275 

.137 

.283 

.110 

.292 
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H?arUy-mt».ntl«d O-owid . 8o*M, 65** 

fr«Ms usad i 70B06 - 70H10, 70B?7 - 70BJ4 
total aroa 1 55.lBO.AlJ km* 
number of crater* t 9 


OVTGT- A'- ‘ ’ 


fllmnrtyr (k») emly.$>y^.Bftt/Mai (*»**) 


8.157 

.181 X 10 

4.456 

.363 

2.232 

.544 

2.108 

.726 

.830 

.109 X 10 

.827 

.127 

.689 

.145 

.551 

.163 

.253 

.181 
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